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f^ ' Abstract 

■ This paper studies systems of particles following independent random walks and subject to 

TT annihilation, binary branching, coalescence, and deaths. In the case without annihilation. 



such systems have been studied in our 2005 paper "Branching-coalescing particle systems" 
The case with annihilation is considerably more difficult, mainly as a consequence of the non- 
V^ I monotonicity of such systems and a more complicated duality. Nevertheless, we show that 

adding annihilation does not significantly change the long-time behavior of the process and in 
fact, systems with annihilation can be obtained by thinning systems without annihilation. 
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1 Results 

1.1 Introduction 

In [AS05) ■ we studied systems of particles that perform independent random walks, branch binarily, 
coalesce, and die. Our motivation came from two directions. On the one hand, we were driven by the 
wish to study a population dynamic model that is more realistic than the usual branching particle 
systems, since the population at a given site cannot grow unboundedly but is instead controlled by 
an extra death term that is quadratic in the number of particles, which can be interpreted as extra 
deaths due to competition. On the other hand, such systems of branching and coalescing particles 
are known to be dual to certain systems of interacting diffusions, modelling gene frequencies in 
spatially structured populations subject to resampling, mutation, and selection |SU86] . In this 
context, the branching-coalecing particles can be interpreted as 'potential ancestors' |KN97) . 

Apart from this duality, which was known, we showed in |AS05j that our particle systems are also 
related to resampling-selection processes by a Poissonization relation. Moreover, we proved that 
systems started with infinitely many particles on each site come down from infinity (a fact that 
had been proved before, with a less explicit bound, in [DDL90^ ) and that systems on quite general 
spatially homogeneous lattices have at most one nontrivial homogeneous invariant law, which, if it 
exists, is the long-time limit law of the process started in any nontrivial homogeneous initial law. 

In the present paper, we generalize all these results to systems where moreover, with some 
positive rate, pairs of particles on the same site annihilate each other, resulting in the disappearance 
of both particles. This my not seem like it should make a big difference with coalescence, where 
only one particle disappears -and indeed our results confirm this- but from the technical point 
of view annihilation has the huge disadvantage of making the system non-monotone, which means 
that many simple comparison arguments are not available. Some pioneering work on non-monotone 
systems can be found in, e.g., |BG85 , Sud90,iDur91| . Despite progress in recent years, non-monotone 
particle systems are still generally less studied and worse understood than monotone ones. 

As in the case without annihilation, our main tool is duality. In fact, it turns out that systems 
with annihilation are dual to the same Markov process (a system of interacting Wright-Fisher 
diffusions) as those without it, but with a different (and more complicated) duality function. As 
a result, we obtain Poissonization and thinning relations which show, among others, that systems 
with annihilation can be obtained from systems without it by independent thinning. We reported 
these duality and thinning relations before (without proof) in |Swa06a) . 

The paper is organized a follows. In Section If .21 we define our model and the dual system 
of interacting diffusions. In Section 11.31 we state our duality result and show how this implies 
Poissonization and thinning relations. Section 11.41 presents our main results, showing that the 
system started with infinitely many particles comes down from infinity and that systems started in 
a spatially homogeneous, nontrivial invariant law converge to a unique homogeneous invariant law. 
Section [1.51 contains more discussion and an overview of our proofs, which are given in Section [21 

1.2 Definition of the models 

Let A be a finite or countably infinite set and let q{i,j) > (i, j e A, i ^ j) be the transition rates 
of a continuous time Markov process on A, the underlying motion, which jumps from site i to site j 
with rate q{i,j)- For notational convenience, we set q{i,i) := (i £ A). We assume that the rates 
q{i,j) are uniformly summable and (in a weak sense) irreducible, and that the counting measure 
on A is an invariant law for the underlying motion, i.e.: 

(i) sup^g(i,j) < oo, 

i 

3 

(n) VA C A, A 7^ 0, A 3i e A, j e A\A such that q{i,j) > or q{j,i) > 0, (i,i) 

(iii) ^qHhj) = ^q{-i,j) Vi e A, where q^{i,j) := q{j,i)- 
j j 

Here and elsewhere sums and suprema over i,j always run over A, unless stated otherwise. 



Branching- annihilating particle systems. We now let (A,g) be as above, fix rates a,b,c,d > 0, and 
consider systems of particles subject to the following dynamics. 

1° Each particle jumps, independently of the others, from site i to site j with rate q(i,j). 

2° Each pair of particles, present on the same site, annihilates with rate 2a, resulting in the 
disappearance of both particles. 

3° Each particle splits with rate b into two new particles, created on the position of the old one. 

4° Each pair of particles, present on the same site, coalesces with rate 2c, resulting in the creation 
of one new particle on the position of the two old ones. 

5° Each particle dies (disappears) with rate d. 

Let Xt{i) denote the number of particles present at site z S A and time t > 0. Then X = (Xt)t>o, 
with Xt = {Xt{i))i£A, is a Markov process with formal generator 

Gf{x) :=^ g(z, j)xW{/(x + S, - S,) - f{x)} + a^x(i)(a;(z) - l){f{x - 2,50 " /(^)} 

ij i 

+ bJ2^imi^ + S^) - fix)} +cJ2<^)«^) - mfix - S^) - fix)} (1.2) 

+dY,xii){fix-6,)-fix)}, 

i 

where Siij) := 1 if i — j and 5iij) := otherwise. We call X the iq, a, b, c, d)-branco-process. 

The process X can be defined for finite initial states and also for some infinite initial states in 
an appropriate Liggett-Spitzcr space. Following [LSSlj . we define 

£:^(A) := {x gN'^ : ||a;||.^ < cx)}, with ||x||.^ := ^7,|a;(i)|, (1.3) 

i 

where 7 = i'Ji)ieA sue strictly positive constants satisfying 

^7, <oo and ^iqiij) + qHiJ))lj < KJ^ (i G A) (1.4) 

i 3 

for some K < cxd. (Our assumptions on q imply the existence of a 7 satisfying (jl.4[) .) 

Resampling selection processes. Let (A, q) be as before, let r, s, m be nonnegative constants, and 
let X = iXt)t>o be the [0, IJ'^-valued Markov process given by the unique pathwise solutions to the 
infinite dimensional stochastic differential equation (SDE) (see |SU86[[AS05) ): 

dXtii) = V qij, i)iXtij) - M^)) dt + sXtii)il - Xti^)) dt ~ mXtii) dt 

3 ^ (1.5) 

+ ^2rXtii)il - Xtii))dBtii) (t > 0, ie A), 

where (i3(i))igA is a collection of independent Brownian motions. The process A" is a system of 
linearly interacting Wright-Fisher diffusions, also known as stepping stone model, which can be 
used to model the spatial distribution of gene frequencies in the presence of resampling, selection, 
and mutation. Following AS05; , we call X the resampling-selection process with underlying motion 
(A,g), resampling rate r, selection rate s, and mutation rate m, or shortly the iq,r, s,m)-resem- 
process. 



1.3 Duality, Poissonization, and thinning 

We start with some notation. For 4't'4' ^ [~oOj oo] , we write 

(</.,V):=^^«VX») and |0| :- ^ |0(z)|, (1.6) 

i i 

whenever the infinite sums are defined. For any </> : A — )■ [—1,1] and a; : A — >■ N we write 

0^:=J|0(i)^« with 0°:=1 (1.7) 

i 

whenever the infinite product converges and the Umit does not depend on the order of the coordi- 
nates. The following proposition generalizes |AS05[ Theorem 1 (a)]. 

Proposition 1 (Duality) Assume that a + c > and let 

a — a/{a + c), r = a + c, s = [1 + a)b, and m = ab + d, (1-8) 

or equivalently 

a — ar, b — s/(l + a), c — (1 — a)r, and d — m ~ as/(l + a). (1-9) 

Let X be a {q, a,b, c, d)-branco-process with Xq g £-y{-^) a.s. and let X^ be a {q\r, s,m)-resem- 
process, independent of X. Suppose that one or more of the following conditions are satisfied: 

(i) a < 1, (ii) |Xol < oo a.s., (in) \X^\ < oo a.s. (1-10) 

Then 

E[(l-(l + a)A-J)^'] =E[(l-(l + a)A-/)^o] (^ > q), (l.H) 

where the infinite products inside the expectation are a.s. well-defined. 

Proposition [U together with a self-duality for {q, r, s, m)-resem-processes described in |AS051 Theo- 
rem 1 (b)] , implies that (q, a, b, c, (i)-branco-processes can be obtained as Poissonizations of resamp- 
ling-selection processes, and as thinnings of each other, as we explain now. (These thinning relations 
will prove useful several times in what will follow. On the other hand, we have no application of 
the Poissonization relations, but since they are very similar and closely related, we treat them here 
as well.) 

If (/) is a [0, oo) -valued random variable, then by definition a Poisson measure with random 
intensity cj) is an N -valued random variable Pois(0) whose law is uniquely determined by 

E[(l-V)^°is('^)] =E[e~('^'^)] (^-£[0,1]^), (1.12) 

where we allow for the case that e~'^'^/ = e~°° := 0. In particular, if (f) is nonrandom, then the 
components (Pois((/i)(i))igA are independent Poisson distributed random variables with intensity 

If X and (p are random variables taking values in N and [0, 1] , respectively, then by definition 
a 0-thinning of x is an N -valued random variable Thin^ (x) whose law is uniquely determined by 

E[(l - ip)^^^Mx)] ^ E[(l - ^^pf] (V- e [0, 1]^). (1.13) 

In particular, when x and (f> are nonrandom and x = ^^ di^ , then a 0-thinning of x can be 
constructed as Thin0(a;) := ^^ Xn^i^ where the Xn are independent {0, l}-valued random variables 
with P[xn = 1] = 4'{in)- More generally, if x and (j) are random, then we may construct Thin0(a;) 
in such a way that its conditional law given x and </> is as in the deterministic case. It is not 
hard to check that (|1.13p holds more generally for any ip £ [0,2]^^ provided (1 — ip'J^'^^^'I'y-'^) is a.s. 
well-defined. 



Proposition 2 (Poissonization and thinning) Fix s,m > 0, r > 0, and < (3 < a < I such 
that m — jTgS > 0. Let X and X be the {q,ar, jr^s, (1 — a)r,rn — j^s)-branco-process and 

(g, /3r, YTflS, (1 — I3)r,m — YTns)-hranco-process, respectively, and let X he the {q,r,s,m)-resem- 
process. Then 

¥[Xo e •] = p[Pois((j^A'o) G •] implies P[X, e •] = p[Pois((^^A',) e •] (t > 0). 

(1.14) 
and 

P[Xo e • ] = P[Thin_i+£ (Xo) e ■ ] implies P[Xt e ■] = P[Tliini+£ (Xt) G • ] (^ > 0). 

(1.15) 

Proof Formula (|1.14l) has been proved in case a = in [AS05| . The general case can be derived 
along the same lines. Alternatively, this can be derived from the case a = using the fact that 
P[Thin_L(Pois(f A:'t)) G ■] = P[Pois( ,^^^^^ A't) e ■], and formula (HIS)) , which we prove now. 

If the initial laws of X and X are related as in (|1.15p and X^ is a (g^,r, s, m)-resem-process 
started in Xq = (f> with |(/)| < oo, then by (|l.lip . 

E[(l - (1 + a)^)^''"^^''*^] = E[(l - (1 + /3)0)^*] - E[(l - (1 + f3)X^f''] 

= E[(l - (1 + a)^'/)^'''"^^''"'] = E[(l - (1 + a)xlf'>] = E[(1 - (1 + a)</.)^'] (t > 0), 

(1.16) 
where we have used that by (ASPS) Lemma 20] one has \Xf\ < cxd a.s. for each t > 0, which 
guarantees that the infinite products are a.s. well-defined. Since (jl.lSp holds for all (p G [0, l]'^ with 
\(j)\ < oo, (fTTSl) follows. I 

As an immediate corollary of formula (|1.15l) . we have: 

Corollary 3 (Thinnings of processes without annihilation) Let a,b,c,d> and a + c > 0. 
Let X be the {q, a, b, c, d)-hranco-process, a :~ ^^, and let X be the {q, 0, (1 + a)b, a + c, ab + d)- 
branco-process. Then 

¥[Xo £■]= P[Thin_L (Xo) e • ] implies P[Xt G • ] = P[Thin_L (Xt) £ ■] (i > 0). (1.17) 

In particular, each branco-process with a positive annihilation rate can be obtained as a thinning 
of a process with zero annihilation rate. 

1.4 Main results 

Let N = N U {cx)} denote the one-point compactification of N, and equip N with the product 
topology. We say that probability measures Vn on N converge weakly to a limit v, denoted as 
Un => ly, when Ji/nidx)f{x) -^ J h'{dx)f{x) for every / e C(N ), the space of continuous real 

functions on N . 

Our first main result shows that it is possible to start a {q, a, b, c, (i)-branco-process with infinitely 
many particles at each site. We call this the (g, a, b, c, (i)-branco process started at infinity. This 
result generalizes |AS051 Theorem 2]. For branching-coalescing particle systems on li'^ with more 
general branching and coalescing mechanisms, but without annihilation, a similar result has been 
proved in |DDL90j . 

Theorem 4 (The mEiximal process) Assume that a + c> Q. Then there exists an £^{A)-valued 
process X^°°' = (Xj )t>o with the following properties: 

(a) For each e > 0, {X^ )t>e is the {q, a, b, c, d)-branco-process starting in X^ ■ 



(b) Set r := a + b + c — d. Then 



ifr^O, 



I, (2a + c)t 

(c) If X^"' are {q,a,b,c,d)-branco-processes starting in initial states a;*-"' G £j{^) such that 

a;(")(i)too as n 1 00 (i £ A), (1.19) 

then 



(d) There exists an invariant measure V of the {q, a,b, c, d)-branco-process such that 



£(X^t">) =^ C{Xr') {t > 0). (1.20) 

n—yoo 



c{x^r^) =^ V. (1.21) 

(e) The measure V is uniquely characterised by 

F(da;)(l - (1 + a)(f)f = P"*[3t > such that x} = 0] (0 G [0, 1]'^, \(j)\ < 00), (1.22) 



where a := a/{a + c) and X'^ denotes the {q\a + c, (1 + a)b, ab + d))-resem-process started in cj). 

(f ) Ifr,s, TO, a, /3 are as in Proposition\^ and X^^^' and X are the corresponding branco-processes 
started at infinity, then 

P[Xt^°°^ G-]=P[Thini+c(x|°°^)e-] (t>0). (1.23) 

A similar thinning relation holds between their long-time limit laws. 

If a = 0, then it has been shown in [AS051 Theorem 2 (e)] that 17 dominates any other invariant 
measure in the stochastic order, hence V can righteously be called the upper invariant measure of 
the process. In the general case, when we have annihilation, we do not know how to compare V 
with other invariant measures in the stochastic order, and we only work with the characterization 
of 77 in ([L221) . 

To formulate our final result, we need some definitions. Let (A, q) be our lattice with jump 
kernel of the underlying motion, as before. By definition, an automorphism of (A, q) is a bijection 
(7 : A ^ A such that q{gi,gj) = q{i,j) for all i,j G A. We denote the group of all automorphisms 
of (A, q) by Aut(A, q). We say that a subgroup G C Aut(A, q) is transitive if for each i, j G A there 
exists a (7 G G such that gi ~ j. We say that (A, g) is homogeneous if Aut(A,(7) is transitive. We 
define shift operators Tg : N^ ^ N^ by 

Tgxij) := x{g-'j) (i G A, x G N^, 5 G Aut(A, q)). (1.24) 

If G is a subgroup of Aut(A, q), then we say that a probability measure i^ on N is G -homogeneous if 
voTg^ = V for all g ^ G. For example, if A = Z'' and q{i,j) = l{|i_j|=i} (nearest-neighbor random 
walk), then the group G of translations i j-t- i + j (j G A) is a transitive subgroup of Aut(A, q) and 
the G-homogeneous probability measures are the translation invariant probability measures. 
The next theorem, which generalizes [AS05[ Theorem 4 (a)], is our main result. 

Theorem 5 (Convergence to the upper invariant measure) Assume that (A, g) is infinite 
and homogeneous, G is a transitive subgroup o/Aut(A, g), anda + c > 0. Let X be the {q,a,b,c,d)- 
branco process started in a G-homogeneous nontrivial initial law C{Xq). Then C{Xt) ^V ast ^ cxd, 
where v is the measure in lll.22\) . 



1.5 Discussion and outline 

The dualities in Proposition[l]and |AS05[ Theorem 1 (b)], as well as the Poissonization and thinning 
relations in Proposition [2] play a central role in the present paper. These relations, whose discovery 
was the starting point of the present work, are similar to duality and thinning relations between 
general nearest-neighbor interacting particle systems discovered by Lloyd and Sudbury in |SL95i 
ISL97[ [SudOOj . In fact, as has been demonstrated in |Swa06al Prop. 6 and Lemma 7] (see also the 
more detailed preprint of the same paper, [SwaOGbl Prop 4.2 and Lemma 4.3]), our relations can (at 
least formally) be obtained as 'local mean field' limits of (a special case of) the relations of Lloyd 
and Sudbury. In |SL97| , Lloyd and Sudbury observed that quite generally, if two interacting particle 
systems have the same dual (whith a special sort of duality relation as described in that article), 
then one is a thinning of the other. This general principle is also responsible for the Poissonization 
and thinning relations of our Proposition [51 

The thinning relation in Corollary |3] is especially noteworthy, since it allows us to compare non- 
monotone systems (which are generally hard to study) with monotone systems. Also, the thinning 
relation (jl.23[) allows us to prove that the unique nontrivial homogeneous invariant measures of 
{q, ar, tt— s, (1 — a)r, m— j^2_s)-branco-processes are monotone in a (w.r.t. to the stochastic order). 
Such sort of comparison results between non-monotone systems are rarely available. In fact, these 
thinning relations suggest that the ergodic behavior of {q,ar, jr^s, (1 — a)r,m — Y^s)-branco- 
processes (with r,s,m fixed but arbitrary a) and the (r, s, TO)-resem process should all be 'basically 
the same'. 

It does not seem straightforward to make this claim rigorous, however. The reason is that 
Poissonization or thinning can only produce certain initial laws. Thus, an ergodic result for 
resampling-selection processes, as has been proved in |SU86j . only implies an ergodic result for 
branching-annihilating particle systems started in initial laws that are Poisson with random inten- 
sity, and likewise, the ergodic result for branching-annihilating particle systems in JAS05) implies 
our Theorem [5] only for special initial laws, that are thinnings of other laws. 

Our main tool for proving the statement for general initial laws is, like in our previous paper, 
duality. In this respect, our methods differ from those in [DDL90J . which are based on entropy 
calculations, but are similar to those used in, for example, |SU86[ IBDD91[ IAS05[ ISS08| . The 
papers |SU86[ IAS05| are particularly close in spirit. The sort of cancellative systems type duality 
that we have to use in the present paper is somewhat harder to work with than the additive systems 
type duality in |SU86|RS05| . Earlier applications of this sort of 'cancellative' duality can be found 
in |BDD91llSS08] . 

The remainder of this paper is devoted to proofs. Proposition [T] and Theorems H] and [5] are 
proved in Sections 12. 5[ 12.61 and 12. 8[ respectively. 

2 Proofs 

2.1 Construction and approximation 

2.2 Finite systems 

We denote the set of finite particle configurations by A/'(A) := {x e N^ : \x\ < oo} and let 

S{J\f{A)) :^ {/: A/'(A) ^R: |/(x)| < K\x\'' + M for some K,M,k> 0} (2.1) 

denote the space of real functions on A/'(A) satisfying a polynomial growth condition. Recall the 
definition of the operator G from (jl.2p . Generalizing |AS051 Prop. 8], we have the following result. 
Below and in what follows, we let P'^ denote the law of the {q, a, b, c, (i)-branco-process started in x 
and we let E'^ denote expectation with respect to P^. 

Proposition 6 (Finite branco-processes) Let X be the {q,a,b,c,d)-branco-process started in a 
finite state x. Then X does not explode. Moreover, with z^^' :~ z{z -|- 1) • • • (^ -|- fe — 1), one has 

E^[|A:t|<'=>] < |x|<'=>e'=''* (fc = 1,2,..., i>0). (2.2) 



For each f G S{Af{A)), one has Gf G S{J\f{A)) and X solves the martingale problem for the 
operator G with domain S{J\f{h)). 

Proof The proof of |AS05i Prop. 8] carries over without a change. I 

We equip N with the componentwise order, i.e., for two states x,x ^ N , we write a; < i if 
x{i) < x{i) for all i G A. In AS0 5 . we made extensive use of monotonicity of branching-coalescing 
particle systems. For systems with annihilation, most of these arguments do no longer work. In 
fact, we can only prove the following fact. 

Lemma 7 (Comparison of branco-processes) Let X and X be the {q,a,b,c,d)-branco-process 
and the {q, 0, 6, c, d)-branco-process started in finite initial states x and x, respectively. Assume that 

x<x, b<b, a + c>c, d > d. (2.3) 

Then X and X can be coupled in such a way that 

Xt <Xt it> 0). (2.4) 

Proof This can be proved in the same way as [AS05[ Lemma 9] , by constructing a bivariate process 
{B, W), say of black and white particles, such that X = B are the black particles and X = B + W 
are the black and white particles together, with dynamics as described there, except that each pair 
of black particles, present at the same site, is replaced with rate 2(1 — 9)c by one black and one 
white particle, with rate 2(1 — 0)a by two white particles, with rate 20c by one black particle, and 
with rate 29a by one white particle, where :~ c/{a + c). I 

We will often need to compare two {q, a, b, c, d)-branco-processes with the same parameters but 
different initial states. A convenient way to do this is to use coupling. Let (Y'^^ ,Y^^ ,Y^'^) be a 
trivariate process, in which particles jump, die and give birth to particles of their own type, and 
pairs of particles of the same type annihilate and coalesce in the usual way of a (g, a, 6, c, ci)-branco- 
processes, and in addition, pairs of particles of different types coalesce to one new particle with a 
type that depends on its parents, according to the following rates: 

01 + 10 h-> 11 at rate r, 

flO at rate 2a, (2.5) 

01 + 11 h^ 111 + , o 

[11 at rate 2c, 

and similarly 10 + 11 iH> 01 or 11 at rate 2a resp. 2c. Then it is easy to see that, for any choice of 
the parameter r > 0, both X := Y^^ + Y^^ and X' := Y^° + Y^^ are {q, a, b, c, (i)-branco-processes. 
We will call this the standard coupling with parameter r. Note that if a = 0, then Xq < Xq implies 
Xt < X'f. for alH > but the same conclusion cannot be drawn if a > because of the transition 
01 + 11h^ 10. 

Let X be the (g, a, b, c, (i)-branco-process. It follows from Proposition [6] that the semigroup 
{St)t>o defined by 

Stfix) := E-[f{Xt)] {t>0, xe Af{A), f G S{Af{A))) (2.6) 

maps S{J\f{A)) into itself. The semigroup gives first moments of functions of our process. We will 
also need a covariance formula for functions of our process, the general form of which is well-known. 
Below, for any measure fi and function /, we write /x/ := J fdfi whenever the integral is well- 
defined, and we let Cov^(/, g) :— ^i{fg) — (m/)(m.'7) denote the covariance of functions /, g under /i. 
Note that if /x is a probability measure on N{A), then ^Stf = / ^(da;)E^(/(Xt)], i.e., ^iSt is the 
law at time t of the (g, a, 6, c, d)-branco-processes started in the initial law /i. 

Lemma 8 (Covariance formula) Let {St)t>Q be the semigroup defined in 12. g|) and let fi be a 
probability measure on M{A) such that j ii{dx)\x\^ < cx) for all fc > 1. Then, for each t > and 
f,g £ S{N{A)), one has 



tg) + 2 / 
Jo 



Coy ^sAf,9) = GoY^iStf, Stg) + 2 / f,St-MSsf, S,g)ds, (2.7) 



where T{f,g) ■- j{G{fg) - {Gf)g - .f{Gg)) is given by 

2r(/, g){x) - Y. QihjH^){f{x + Sj - S,) - fix)) {g{x + S, - S,) - g{x)) 

ij 

+a^a;(i)(a;(i) - l)(/(x - 26,) - f{x)){g{x - 26,) - g{x)) 

i 

+bJ2xii){f{^ + Si) - fix)) {g{x + 6i) - gix)) (2.3) 

+c^a;(z)(a;(z) ~ l){f{x - 6,) - fix)) {g{x - 6,) - g{x)) 

i 

+d^x(i)(/(x - 6,) - fix)) igix - 6,) - gix)). 

i 

Proof Formula (|2.7p is standard, but the details of the proof vary depending on the Markov 
process under consideration. In the present case, we can copy the proof of [Swa09[ Prop. 2.2] 
almost without a change. We start by noting that fg G 5(A/'(A)) for all f,g G 5(A/'(A)), hence 
Tif,g) := i(G(/5) - (G/)g - fiGg)) is well-defined for all f,ge 5(A/'(A)). It is a straightforward 
excercise to check that r(/, g) is given by (|2.8p . Now (|2.7p will follow from a standard argument 
(such as given in jSwa09[ Prop. 2.2]) provided we show that 

^^Ss{iStf)iSug))^S,G{iStf)iSug)), 

|5,((5J)(^„5))=^.((G5J)(^n.9)), (2.9) 

£S,{iSJ)iS,,g))=Ss{iSJ)iGS,,g)) 

for all < s,t,u and f,g G 5(7V(A)). Let us say that a sequence of functions /„ e 5(A/'(A)) 
converges 'nicely' to a limit / £ 5(A/'(A)) if /„ — > / pointwise and there exist constants K, Ad, k > 
such that sup„ |/n(a;)| < iir|a;|'^' + M. Then (|2.2p and dominated convergence show that fn —^ f 
'nicely' implies Stfn -> Stf 'nicely'. Note also that if /„,/, 5 G 5(A/'(A)) and /„ ^- / 'nicely', 
then fng -^ fg 'nicely'. It is easy to check that Gf G 5(7V(A)) for all / G 5(7V(A)). Since the 
(g, a, &, c, d)-branco-process X^ started in a deterministic initial state X§ = a; G A/'(A) solves the 
martingale problem for G, we have 



r'{Stfix)-fix))=t-' f E[Gfix:)]d.s^Gfix) (xGA/-(A)), (2. 

Jo *^'-' 
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which by (g^) and the fact that Gf G 5(A/'(A)) implies that t^^iStf - f) ^ Gf 'nicely' as i 4, 0. 
Combining three facts, we see that 



§-^Ss{iStf)iSug)) = ^^SsiPe - l)iiStf)iSug)) = 5,G((5J)(^„5)), 
f 5,((5,/)(5„<?)) = hm5,(((F, - l)5t/)(5„5)) = 5,((G5J)(5„g)), 



and similarly for the derivative w.r.t. u, where we are using that if the right-hand derivative of a 
continuous real function exists in each point and depends continuously on i, then the function is 
continuously difFerentiable (see, e.g., |HS75[ Excersise 17.24]). I 

2.3 Infinite systems 

Recall the definition of the Liggett-Spitzer space £^7 (A) from (|1.3p . We let Ci^ipiS^iK)) denote the 
class of Lipschitz functions on £7 (A), i.e., / : E-^i^) -^ K such that |/(x) — fiy)\ < L\\x — y\\.y for 
some L < 00. 

The main result of this section is the following generalization of [AS05| Prop. 11]. 

Proposition 9 (Construction of branco-processes) Let iSt)t>o be the semigroup defined in 
^2.6\) . For each f G Ci^i-piE^iK)) and t > 0, the function Stf defined in i2.6\) can be extended to 



a unique Lip schitz Junction on £~f{A), also denoted by Stf ■ There exists a unique (in distribution) 
time-homogeneous Markov process with cadlag sample paths in the space ^'y(A) equipped with the 
norm \\ ■ ||-y, such that 

E-[/(X,)] = Stf{x) (/ e CLip(f^(A)), X e £^{A), t > 0). (2.12) 

To prepare for the proof of Proposition [9l we start with the following lemma, which generalizes 
(lEUSl Lemma 12]. 

Lemma 10 (Action of the semigroup on Lipschitz functions) Let {St)t>o be the semigroup 
of the {q, a,b, c, d)-branco-process, defined in i2.6\) . If f : A/'(A) -^ S. is Lipschitz continuous in the 
norm \\ ■ \\-y from ^1-4^ , with Lipschitz constant L, then 



\Stf{x)-Stf{x')\<Le^''+'-''^'\\x-x'\\^ {x,x' eMiA), t>0), (2.13) 

where K is the constant from \1.4^ . 

Proof Let X ~ Y^^ + F^^ and X' — Y^'^ + Y^^ be (q, a, b, c, (i)-branco-processes started in Xq = x 
and Xq = x' , coupled using the standard coupling from (|2.5p . in such a way that {Yq^ ,Yq^ ,Yq'^) = 
{{x — x')+,x A x' , {x' — x)^). Then 

\Stfix) - Stf{x')\ - \nfiXt)] - E[fiX',)]\ < E[\fiXt) - fiX^)\\ 
<LE[\\X,-X^\\,]=LE[\\Y°' + Y,'%]. 

Let us choose the parameter r in the standard coupling as r :— 2{a + c). Then it is easy to see 
that (F°^,y^°) can be coupled to a {q,0,b,a + c, d)-branco-process Z started in Zq — \x — x'\ in 
such a way that Y^°'^ + Y^^^ < Zt for all t > 0. Therefore, by |AS051 formula (3.13)], we can further 
estimate the quantity in (|2.14l) as 

\Stfix) - Stf{x')\ < LE[\\Zt\\^] < Le'^^+'-^^x - x'\\^. (2.15) 



Proof of Proposition [9] Since J\f{A) is a dense subset of ■^^(A), Lemma [TOl implies that for each 
/ G Cup{£-i{A)) and t > 0, the function Stf defined in (j2.6p can be extended to a unique Lipschitz 
function on £j{A). The proof of Lemma [TOl moreover shows that two {q, a, b, c, (i)-branco- processes 
X, X' started in finite initial states x, x' can be coupled such that 

E[||Xt-X;||^] <e(^'+''-'^)*||a;-a;'||^ (i > 0). (2.16) 

It is not hard to see that for each x G ^7 (A) we can choose Xn G A/'(A) such that IIj:^ — ^|| ^ and 

00 

^\\Xn~ Xn-l\\^ < 00. (2-17) 

n=l 

(For example, any x„ t 2; has these properties.) Let X" be the process started in Xq — x„. By 
(|2.16p . we can inductively couple the processes X°, X^,X^, ... in such a way that 

E[\\X^ - Xr%] < e(^+''-'')*||x„ - Xn-ih {n>l,t> 0). (2.18) 

It follows that for each (deterministic) i > 0, the sequence (X")"-*' is a.s. a Cauchy sequence in 
the complete metric space £j{A), hence for each t > there a.s. exists an £^..y (A)-valued random 
variable Xt such that jjAT" — Xt\\^ — > 0. By Fatou, 

00 
E[||Xr - Xt\U] < liminf E[||Xr - Xr\\^] < e^^+b-d)* V ||x„+i - x„||^ -^ 0. (2.19) 



71— >-oo 
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Just as in |AS05[ Lemma 13], it is now straightforward to check that {Xt)t>o is a Markov process 
with semigroup {St)t>o- Note, however, that in the arguments so far we have only constructed 
X = {Xt)t>Q at deterministic times. To show that X has a version with cadlag sample paths 
(where the limits from the left and right are defined w.r.t. the norm || • ||^), we adapt arguments 
from the proof of |AS05[ Prop. 11]. It suffices to prove X has cadlag sample paths on the time 
interval [0, 1]. 

Let V be the process with generator 

Gvfix) := Y. q{i,j)x{i){f{x + 5,)~ f(x)} + h ^ x(z){/(x + S,) - /(x)}, (2.20) 

ij i 

which describes a branching process in which particles don't move or die, and each particle at i 
gives with rate q(i,j) birth to a particle at j and with rate 6 to a particle at i. We claim that 
a (g, a, 6, c, d)-branco-process X, started in a finite initial state Xq — x, can be coupled to the 
process V started in Vq = x in such a way that Xt < Vt for all i > 0. To see this, let {B, W) be 
a bivariate process, say of black and white particles, started in {Bo,Wo) = (a;,0), such that the 
black particles evolve as a (g, a, &, c, ci)-branco-process, the white particles evolve according to the 
generator in (|2.20l) . and each time a black particle disappears from a site i due to jumps, annihilation 
or coalescence, a white particle is created at i. Then it is easy to see that X — B and V = B + W. 
By [AS05. formula (3.25)], 

IE[||V^t||,]<e(^+'')*||x||^, (2.21) 

where K is the constant from (jl.4p . Since V is nondecreasing in t, since Vt{i) increases by one each 
time Xt{i) does, and since X cannot become negative, it follows that 

\{te [0,1] : Xt-ii) ^ Xtii)}\ <xii)+2Viii). (2.22) 

Applying this to the process X", multiplying with 7^ and summing over i, we see that 

5]7.IE[|{ie[0,l]:Xr_(*)y^Xr(z)}|] <(l + 2e^'+^)||a;„||^, (2.23) 

i 

which by the convergence of l|a;„||-y gives us a uniform bound on the number of jumps made by X"'. 
We wish to show that for large n, the processes X" and X"+^ make mostly the same jumps. 
To this aim, let X" = Y°^ + F" and X"+i = Y'^° + y" be two (g,a, &, c,d)-branco-processes, 
coupled using the standard coupling from (12.51) . with r = 2(a + c) and {Yq^,Yq^,Yq^) = ((a;„ — 
Xn+i)+,Xn A Xn+1, {xn+1 ~ 2^)+)- Thcu, just as in the proof of Lemma [TUl the process (y°-^,y^°) 
can be coupled to a (g, 0,6,a + c, (i)-branco-process Z started in Zq = |x„ — a;„+i| in such a way 
that Y^^ + Y^'^ < Zt for all t > 0. Likewise, it is not hard to see that we can couple {Y'^^, Y^'^) to a 
process V with dynamics as in (|2.20l) started in Vq = |a;„ — a;„+i | , in such a way that Y^^ + Y^^ < Zt 
for all i > and moreover, whenever Y^^{i) or Y'^^{i) increases, the process V{i) increases by the 
same amount. Let 

J„(z) ■.= {te [0, 1] : XlUt) ^ X^t)} (2.24) 

be the set of jump times up to time one of the process X"(i) and let 

/(z) := {t E [0, 1] : r,Oi(z) ^ Y^^^i)} U {t e [0, 1] : r,"'(*) ^ y,i°(*)}. (2.25) 

Then the symmetric difference J„(i) A Jn+i{i) — (>/n(i)\>^n+i(i)) U {Jn+i{i)\Jn{'i)) of J„(i) and 
Jn+i{i) is contained in I{i) and, by the arguments leading up to (|2.22p . \I{i)\ < \xn{i) — Xn+i{i)\ + 
2^1(1). Thus, in analogy with (|2.23p . we find that 

^7iE[|J„(i)AJ„+i(i)|] <(l + 2e^+^)|lx„-x„+i||^. (2.26) 

i 

By (|2.17p . it follows that the sets Jn{i) converge as n ^- c», i.e., for each i G A there is a (random) 
n such that J„(i) = Jn+i{i) = Jn+2{i) = ■ ■ ■ ■ Taking into account also (|2.22l) . it follows that the 
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limit process {X{i))t>o has cadlag sample paths for each i e A and the set of jump times of X^{i) 
converges to the set of jump times oi X{i). The fact that the sample path of {X)t>o are also cadlag 
in the norm || • ||^ can be proved in the same way as [AS05i formula (3.31)]. I 

The proof of Proposition [S] yields a useful side result. 

Corollary 11 (Approximation with finite systems) Let x £ ^-yi-^) and a;„ e -^(A) satisfy 
\\xn ~x\\-y — !■ and J2n>i ll^n ~ Xn-1 \\-f < cxD. Then the {q, a, b, c, d)-branco-processes X", X started 
in Xq = Xn and Xq = x can be coupled in such a way that ||X" — Xt\\j -^ a.s. for each t > 0. 

2.4 Covariance estimates 

In this section, we give an upper estimate on the covariance of two functions of a {q,a,b,c,d)- 
branco-process, which shows in particular that events that are sufficiently far apart are almost 
independent. 

For any continuous / : £7 (A) -^ M, we define Sf : A ^ [0, cxd] by 

df{i):^ sup \f{x + S,)-f{x)\ (leA). (2.27) 

It is easy to see that for each continuous / : 6^ (A) — J> R, 

|/(x)-/(2/)|<^5/(z)|x(*)-y(*)| {x,yeS^{A)). (2.28) 

i 

Lemma 12 (Lipschitz functions) A continuous function f : S-yiA) — >■ M is Lipschitz with respect 
to the norm \\ ■ \\j if and only if there exists a constant L < 00 such that Sf{i) < L^i (i £ A). 

Proof If / G CLip(^7(A)), we have \f{x + Si) — f{x)\ < L\\{x + Si) — 2;||^ = L'fi, where L is the 
Lipschitz constant of /, hence Sf{i) < Lji {i e A). Conversely, if the latter condition holds, then 
by (12251) 

\fix)-fiy)\<Lj2iM^)~yi^)\^L\\^~yh (x,j/ef^(A)). (2.29) 



We let B~f (A) denote the space of all functions (/) : A ^> M such that 

sup7,"V(*)l < oo- (2-30) 

i 

Note that by Lemma [H Sf G 5^ (A) for each / G Cupi^-yiA)). 

Let Pt{i,j) denote the probability that the random walk on A that jumps from k to / with rate 
q{k, I), started in i, is a time t located at the position j. For any cj) G By{A), we write 

Ptcb{i):^Y.Pt{^,j)m (i>0, zgA). (2.31) 

j 

It is not hard to check that Pt is well-defined on B.y{A) and maps this space into itself. 

Recall that {St)t>o denotes the semigroup of the (q, a, 6, c, (i)-branco-process, defined in (|2.6p . 

Lemma 13 (Variation estimate) For any {q,a,b,c,d)-branco-process, one has 

SStf < e^'-'^^'PtSf (i > 0, / G Cup{£^{A))) . (2.32) 

Proof Fix i G A and let X = Y°^ + F" and X' = Y'^° + F" be (g, a, 6, c, d)-branco-processes 
started in Xq ~ x and ATq = a; + (5^, coupled using the standard coupling from (j2.5p . in such a way 
that {Y§\Y,}\Y^'^) = {0,x,S^). Then 

\Stf{x) - St fix + S,)\ - \E[fiXt)] - E[/(X;)]| < E[\fiXt) - fiX',)\] 

<E[J2Sfij)\Xtij)~Xiij)\]^J2^f^^My'\j) + Yt'"ij)]<J2^fij)''^'''^'P^('^j)^ 

3 J J 

(2.33) 
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where in the last step we have used that Y'^^ + Y^'~' can be estimated from above by a {q, 0, b, 0, d)- 
branco-process. I 

Proposition 14 (Covariance estimate) Let X — (Xt)t>o be a {q,a,b,c,d)-branco-processes 
started in Xq = x G f-y(A). Then, for each t > 0, there exist functions Kf : A x A^ ^ [0,oo) 
and Lt : A^ X A^ ^ [0, oo) satisfying 



an 



d sup 2, ^ t{i'T k ^ Q) < oo^ and sup >^ Lt(i, j; fc,0) < oo (T < oo), (2.35) 



such that 



\Goy,{f{Xt), g{Xt)) \ < ^ x{z)Kt{i; k, l)Sf{k)Sg{l) 

i,k,l 

+ J2 x{i)xU)Ltit,j;k,l)Sfik)Sgil). 

i,j,k.l 



(2.36) 



for all bounded functions f,g ^ Ci^ip{£.y{A))^ . 



Proof It suffices to prove the claim for finite initial states x e 7V(A). For once the proposition is 
proved for finite systems, for arbitrary x G ^^(A) we can find Af{A) 3 Xn t x- Then by CorollarvfTTl 
the processes X'^,X started in Xn,x can be coupled such that ||Ar" — X^H^ — > for each i > 0, 
hence by bounded pointwise convergence, the left-hand side of (j2.36p for X" converges to the same 
formula for X, while the right-hand side is obviously continuous under monotone limits. 

We will show that for finite systems, the estimate ()2.36p holds even without the boundednes 
assumption on f,g. We apply Lemma [5] A little calculation based on (|2.8p shows that 

2|r(/, g)ix)\ < J2 q{i.j)x{i){5f(i) + 6fij)) {Sg{i) + 6gij)) 

+ {2a + c)J2 x{i){x{i) - l)5f{^5g{i) (2.37) 

■i 

Hh + d)Y,x{i)5f(i)5g(i). 

i 

In view of Lemma [T3l we define P* := e^''~'^^*Pt. Then ^1^, ([237)) and Lemma fT3l show that for 
processes started in a deterministic initial state, 

\CoY4f{Xt),g{Xt))\ 
< I J2'i('^^^iP^^f(')+P^^f'^j))iP^^9{i)+PsSg{j))E-[Xt-s{ms 

■Jo tj 

/■* (2 38) 

+ i2a + c) y2PsSfii)Pjgii)E-[Xt-sii)iXt^sii)-l)]ds 

t ' 

+ {b + d) f y2Ps6f{t)Ps6g{i)E''[Xt-s{i)]ds, 
Jo , 

Let Y = {Yt)t>o be the (g, 0, 6, 0, (i)-branco-process started in Yq = x. By Lemma[71 we can couple 
X and Y such that Xt < Yt for alH > 0. We estimate 

(i) E^[Xti^)]<E^[Yt{^)] = ^ a;(j)A(j, *), 

j (2-39) 

(ii) E-[Xt{i){Xt{i) - l)]<E-[Yt{if] ^W[Ytit)]^ +Yav.AYt{i)). 
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To estimate Yaix{Yt{i)), we apply (|2.38p to the process Y and f = g '■= fi where fi{x) :— x{i). 
Since the annihilation and coalescence rates of Y are zero, this yields 



" ijkl 



+ (2a + c) / ^ x{i)xiJ)Pt-s(.hm)Pt.s{j,m)B,{m;k,l)Sfik)Sgil)ds 

ijklm 

+ (2a + c)/dsV / dux{i)Pu{iJ)CuU;m,m)B,{m;k,l)Sf{k)Sg{l). 



(2.40) 



(2.41) 



Var,(yt(z)) 

< / E 9(^' k){PsSMj) + PsSMk)) (PsSMj) + PsSMk))E-[Yt-s{j)]ds 

+ {b + d) I y2Ps6Mj)Ps6Mj)E^[Yt-s{j)]ds. 
Jo ^ 

Define 

^(z;fc,Z):=E<z(*,j)(A(*,fc) + P,(j,fc))(P,(*,0+Pt(j,0), 

j 
Bt{i;k,l):=Pt{i,k)Pt{iJ). 

Then (|2.38l) can be rewritten as 

\Cov.4f{Xt),g{Xt))\ < f Y.E-[Xt-s{i)]{As{i;k,l) + {b + d)Bs{i;k,l))6f{k)Sg(l)ds 

+ {2a + c) f Y,E^[Xt^s{t){Xt-s{i) - l)]Bs{r,k,l)Sf{k)Sg{l)ds, 
•^0 ^kl 

while (|2.40l) can be rewritten as 

Var,(r,(*))< / EE-[r,_,(j)](^(j;z,z) + (6 + d)i?,(j;*,z))ds, 
Jo ^ 

where we have used that Sfi{j) — l{i=j}- Setting 

Ct{i; fc, := A,(z; fc, /) + (& + d)P,(z; fc, 0, (2.44) 

and inserting (|2.39p and (|2.43l) into (|2.42l) . we obtain 

ICoY^f {Xt),g{Xt))\ 

< f Y.x{t)Pt-s{i,j)CsiJ;k,l)Sfik)dgil)ds 



(2.42) 



(2.43) 



(2.45) 



ijkl 771 

Recalling the definition of Bt{i]j^ k)^ this shows that f|2.36p is satisfied with 

Ktii;k,l):^ f dsy2Pt-,{t,j)Cs{j;k,l) 
Jo ^ 

pt pt — S 

+ {2a + c) / dsV / duPu{i,j)Cu{j;m,m)Ps{m,k)Ps{m,l), (2.46) 

Lt{i,j;k,l):={2a + c) [ dsy2Pt^,{i,m)Pt-s{j,m)Ps{m,k)Ps{m,l). 

The invariance of Kt and Lt under automorphisms of (A, q) is obvious from the analogue property 
of Pt, but the summability condition (12.35^ needs proof. Since P((i, • ) is a probability distribution 
and since the counting measure on A is an invariant law for Pt by assumption (jl.ll) (iii), 

EPt(z,j) = l = E^*(^'*) (i>0, iGA). (2.47) 
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Setting \q\ := J^j <lihj) = Ej 90'>«)> we see that 

^A,(z;fc,0 = ^g(*,j)(A(*,fc)+A(j,fc))(A(*,0 + A(j,/)) 

= ^g(^,J)2e(^-'^)*(A(^,^) + A(J■0)-4|q|e2(h-'^)* (/ S A), 

and therefore, by a similar calculation for Bt{i] k, I), 

Y,Ct{r,m,m)<Y,Ct(j-k,l)^{A\q\+h + d)e'('-''^' {l,meA), (2.49) 

J jk 

which by (|2.46p implies that 

V Kt{i; k, I) < (4|g| + b + d) f ds e(''-'')(*--)e2(&-'^)- 

+ (2a + c) (4|g| +b + d) f As f du e(*'~'^)"e2(''-'^)"e2(''-'')^ < oo, (2.50) 



^2(b-d){t-s)2(b-d)s 



^Lt{i,j;k,l)<{2a + c) f dse^'^''-''''^'-'>e'^'-''^' <oo (t > 0, I e A). 



Corollary 15 (Exponential functions) Let X — {Xt)t>o be a {q,a,b,c,d)-branco-processes 
started in Xq — x E E-y{K), and let fi : A ^ [0, cxd) satisfy J2i m(*) < oo. Then 

i:^[e~T.^^^i^)^t{i^ _TTE^[e"^(*)^*(*)]| 

< i J2 2^(*)^t(^; k, l)fi{k)fi{l) + \ J2 x{z)x{])LS,r, k, l)Kk)fi{l), (2.51) 

i,k,l i^j,k,l 

where Kt,Lt are as in Proposition \14\ 

Proof We first prove the statement if ^ is finitely supported. Let support(/x) = {fci, . . . , km} and 

set 



/„(a;) := e-^(^"M^") and gpix) :^ Y[ Ux). (2.52) 

Q = l 

Then 

E"[g„,(XO] =E"[g„,_i(XO]E"[/„(XO] +Cov4g™_i(X0,/™(X0) 

= E"[g™_2(^t)]E"[/„_i(Xt)] +Coy4grn-2{Xt)Jrn-i{Xt)) 

+ CoVr,[gm-i{Xt), f„i{Xt)) 



(2.53) 



l[E[UXt)] + J2Cov4gc.-i{Xt)Jc.iXt)). 



a=l ai=2 

Therefore, since 

'^^"^'"^"i otherwise, dU(k) - < q otherwise, ^^'^^^ 
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Proposition [in tells us that 

\E^[g,4Xt)]~l[E[f4Xt)]\ 

rn a — 1 ^ (z.OOJ 

< ^ ^ y^x{i)Kt{i]kp,ka)n{ki3)ii{ka) +^x{i)x{j)Lt{i,j\kii,ka)n{kp)iJ.{ka 

a=2 13 = 1 i i,j 

To generalize the statement to the case that J^i ^(0 < oo but n is not finitely supported, it suffices 
to choose finitely supported /i„ t M ^'^d to observe that all terms in (|2.5ip are continuous in /i w.r.t. 
increasing limits. I 

2.5 Duality and subduality 

Recall the definition of 0^ from ()1.7p . 

Lemma 16 (Infinite products) Let < a < 1, cj) e [0, l]'^, x e N^. 

(a) Assume that one or more of the following conditions are satisfied: 

(i) a < 1, (ii) l^l < oo, (iii) 101 < oo. (2.56) 

Then (l — (1 + q;)0) is well-defined. 

(b) Assume that cj) is supported on a finite set and Xn G N converge pointwise to x. Then 
(l — (1 + a)(f>y" — S- (l - (1 + a)(/))'^ as n ^ oo. 

(c) Assume that \(j)\ < cxd and N 3 Xn t ^- Then (l — (1 + ajcf) " — ;• (l — (1 + a)(j>) as n ^>^ oo. 

(d) Assume that either a < 1 or \(j)\ < oo, and let [0, 1] 3 0„ t 0- Then (l — (1 + a)0„) -^ 
(l — (1 + a)(j)) as n ^ oo. 

Proof Since (l - (1 + «)</))"' := Ui (l - (1 + «)</-(«)) "'^*\ where -1 < (l - (1 + a)(/i(i))"'^'' < 1, 

the only way in which the infinite product can be ill-defined is that Ilj |l ^ (1 + Q;)0(*)| > 

while (l — (1 + a)(j){i)^ < for infinitely many i. If a < 1, then —1 < —a < 1 — (1 + a)(j){i), 

so if (l — (1 + a)(l){i)) < for infinitely many i, then Jli |l ^ (1 + Q;)0(*)| = ^ind the 

infinite product is always well-defined. If \x\ < oo, then (l — (1 + a)(j){i)) = 1 for all but finitely 
many i, hence the infinite product is certainly well-defined. If \(f>\ < oo, finally, then <j){i) > ^ for 

finitely many i, hence (l — (1 -I- a)(/)(i)) < for finitely many i and the infinite product is again 
well-defined. This completes the proof of part (a). 

Part (b) is trivial since all but finitely many factors in the infinite product defining (l — (l+a)(/)) 
are one if cj) is finitely supported. 

To prove part (c), we split the product Yl- (l ^ (1 + '^)4'{'i)) iii finitely many factors where 
4>{i) > i and the remaining factors where 4>{i) < ^ and hence (l — (1 + a)(f){i)^ > 0. Then the finite 
part of the product converges as in part (b) while the infinite part converges in a monotone way. 

For the proof of part (d) set I :— {i G A : x{i) ^ 0} and let /_,/o,/+ be the subsets of 
I where 1 — (1 -f a)(/){i) < 0, = and > 0, respectively. If /q 7^ then it is easy to see that 
(1 — (1 -I- a)(j)n) — > = (1 — (1 + a)(j)) , so from now on we may assume that /q = 0. Note 
that 1 — (1 -I- a)(/)n(i) > 1 — (1 + a)<j>{i) > for all i G I+. Therefore, if /_ is finite, as must 
be the case when \(j)\ < 00, then Yiiei (^ ^ (^ + ct)</'n(*)) converges since /_ is finite while 
11,6/+ (1 - (1 + a)Mi)V'^'^ i Uiei+ (1 - (1 + a)(l3ii)y'''\ If /_ is infinite and a < 1, then 
the fact that |l — (1 + Q;)0„(i)| -^ ((1 + a)(/)(i) ~ l) < a for each i <E I- implies that 
[l - {1 + a)<l),,Y ^ {) = {\ - {1 + a)<t>Y . I 

We equip the space [0, 1]^ with the product topology and let C([0, 1]^) denote the space of continuous 
real functions on [0, 1] , equipped with the supremum norm. By Cgjj([0, 1] ) we denote the space 
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of C^ functions on [0,1]'^ depending on finitely many coordinates. By definition, Cgmjj([0, 1]^) 
is the space of continuous functions / on [0, 1]^ such that the partial derivatives -Q§n\f{4>) and 

dcb(i)d6( •) /('?'') sxist for each x G (0, 1)^ and such that the huictions 

'/'^(4l)/(^)).eA ^"d 0^(^^^^/(0))^^^.^^ (2.57) 

can be extended to continuous functions from [0, 1]^ into the spaces ^^(A) and fi{A'^) of absolutely 
summable sequences on A and A^, respectively, equipped with the £^-norm. Define an operator 
^:Ci^([0,in^C([0,l]A)by 

GfW :-E 9(J' *)('^(^') - <^(*)) Wy/(<^) + ^ E '^(*)(1 - '^(*))4l) /('^) 

y ■' (2 58) 

+r5]^«(l-0(z))^/(<^)-m^0(z)4^/(^) (^£[0,1]^). 

i i 

One can check that for / e Cg^^([0, 1]'^), the infinite sums in (I2.58P converge in the supremumnorm 
and the result does not depend on the summation order |Swa99[ Lemma 3.4.4]. It has been shown 
in |AS05[ Section 3.4] that solutions to the SDE (|1.5p solve the martingale problem for the operator 
Q. In view of this, we loosely refer to Q as the generator of the (g, r, s, TO)-resem-process. 

Proof of Proposition [1] Since by Proposition IH] (resp |AS05[ Lemma 20]), \Xq\ < oo (resp. 
\Xo\ < oo) implies \Xt\ < oo (resp. \Xt\ < oo) for all t > 0, by Lemma [TCI each of the conditions 
p.lOP (i)-(iii) guarantees that both sides of equation (|l.lll) are well-defined. 

It suffices to prove (|1.11[) for deterministic initial states, i.e., we want to prove that either a < 1, 
|a:| < oo, or \(f)\ < oo imply that 

E^[{l~il + a)^)^t]^E^[{l-{l + a)X^f] (t > 0), (2.59) 

where E^ and E''^ denote expectation w.r.t. the law of the process X started in Xq = x and the 
process X started in Xq = (j), respectively. We start by proving (|2.59p if \x\ < oo. We wish to apply 
[AS05| Thm 7]. Unfortunately, the original formulation of this theorem contains an error, so we 
have to use the corrected version in |AS09b|, Corollary 2] (see also [AS09a|, Corollary 2]). We apply 
this to the duality function 

*(x,0):= (1-(1 + «)(/.)" (a;GAA(A), (/)€ [0,1]^). (2.60) 

Since |^(a;, 0)| < 1, we obviously have 'i'{-,(j)) £ S{N'{Aj) for each (f) G [0,1]^. Since for each 
X e A/'(A), the function ^(x, •) depends only on finitely many coordinates, we moreover have 
^(x, ■) G Csu^([0, 1]^) for each such x. Let Gbe the generator of the (g, a, 6, c, (i)-branco-processand 
let t/t denote the generator of the (g^, r, s, ?Ti)-resem-process. In order to apply AS09b| Corollary 2], 
we need to check that 

<^{x, 0) := G*( • , <l>){x) =. g^^{x, ■ )((/.) (x G Af{A), ^ G [0, 1]^) (2.61) 

and 

[ ds f dtE[\^{X,,Xt)\] <oo (r>0). (2.62) 

JQ Jo 

To check (|2.6ip . we calculate: 

G*( ■,cj,){x)^J2 «(*' J>W(1 - (1 + ^)^f' ((1 - (1 + ")^y' - (1 - (1 + «)'^)'0 
+a^a;(z)(a;(z) - 1)(1 - {1 + a)<l>r-''^{l - (1 - (1 + a^^) 
+bJ2^{i){l - (1 + a)<l>r-'^ ((1 - (1 + a)0)2*- -(!-(! + «)</.)*') 



-c V x(z)(x(z) - 1)(1 - (1 + aW-^'^ ((1 - (1 + aW' - (1 - (1 + a^' 



+dJ2^im - (1 + «)</-)""'• (1 - (1 - (1 + c^W')- 



(2.63) 
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Since 



and 



^(1 - (1 + a)0f = -(1 + a)a:W(l - (1 + a)0r-^% 
9%(1 _ (1 + a)c^y = (1 + afx{i){x{i) - 1)(1 - (1 + a)0)"-2^' 






+Y^(i + ")0«(i - (1 + «)</-«) 4iy(i - (1 + «)0r 

r(l + a)</'W(l - (1 + a)^iz))gg^{l - (1 + a)(/.f 



Reordering terms gives 

GvI/( •, </.)(a;) = ^ g(*,j)(0(jV '^(*))4l)(l - (1 + «)'/') 
'2a + c 



-C 



+ ((fe - d)^{i) - 6(1 + a)0(j)2) ^(1 - (1 + «)0)- 

+(« + c)4>{t){i ~ ^{^))gg^{i - (1 + «)0r 

+ (1 + a)6(/)W(l - m}ami^ - (1 + ")'^)' 

-(a6 + d)<^(z)4^(l - (1 + a)0)- - etvi;(^, . )(0), 



where we have used (jl.Sp . which imphes in particular that 

2a + c_ 2a + c _(2a + c)(a + c) 



= a 



(2.64) 



(l-(l + a)(/))^'=l-(l + a)0(z), 
(1 - (1 + a)0)2^'- - (1 - (1 + a)0(*))', (2.65) 

(1 - (1 + a)cb)^' - (1 - (1 + a)</.)2*- = (1 + a)</.(i) (1 - (1 + a)(/.(i)) , 
we can rewrite the expression in p.63p as 

G*( • ,0)(a:) =^g(i,j)(0(jv 0«)a|i)(i - (1 +«)0r 



(2.66) 



(2.67) 



(2.68) 



1 + a l + a/(a + c) (a + c) + a 

It is easy to see from (12. 63^ that there exists a constant K such that 

mx, 0)1 < K{1 + |x|2) (0 e [0, 1]^, a; G Af{A)), (2.69) 

hence (|2.62p follows from Proposition [6l This completes the proof of (|2.59p in case |a::| < oo. 

We next claim that ()2.59p holds if a; G ^7 (A) and <j) is supported on a finite set. Choose 
A/'(A) 5 Xn t X and let X" denote the (9, a, 6, c, d)-branco-process started in Xq = a;„. Then 
Corollary [TT] implies that the X" can be coupled such that X"(«) -^ Xt{i) a.s. for each i G A. 
Therefore, taking the limit in (|2.59l) . using the fact that the integrands on the left- and right-hand 
sides converge in a bounded pointwise way by Lemma [T6l (b) and (c), respectively, our claim follows. 

To see that (|2.59p holds more generally if a < 1 or \(j>\ < cx), we choose finitely supported (j)nt4' 
and let A"" denote the ((7,r, s,?7i)-resem-process started in Xq — (j>n. Then |AS05| Lemma 22] 
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implies that the A"" can be coupled such that XI'' {i) t Xt{i) a.s. for each i £ A. The statement then 
follows by letting n — )> cxd and applying Lemma [T51 fd). I 

Fix constants /3 e M, 7 > 0. Let A^(A) := {0 G [0, 00) : \(j)\ < 00} be the space of finite measures 
on A, equipped with the topology of weak convergence, and let y be the Markov process in A^(A) 
given by the unique pathwise solutions to the SDE 



dyS) = V a(j, i)iyt{j) - ytm dt + pytii) dt + ^2^ytii)dBt{i) (2.70) 



j 



(t > 0, i e A). Then y is the well-known super random walk with underlying motion a, growth 
parameter (3 and activity 7. One has |Daw93[ Section 4.2] 

]E0[e-(3^t,V')] ^ e-(0,WtV') (2.71) 

for any G A^(A) and bounded nonnegative t/; : A — > R, where Ut = Utip solves the semilinear 
Cauchy problem 

f ut(z) = Y. "(^■' ^)MJ) ~ Mi)) + /3wt(i) " 7Wt(i)' (i e A, i > 0) (2.72) 

with initial condition uq — tp. The semigroup (Ut)t>o acting on bounded nonnegative functions ijj 
on A is called the log-Laplace semigroup of y. 

It has been shown in [AS05, Prop. 23] that the (g, a, 6, c, (i)-branco-process and the super ran- 
dom walk with underlying motion q^ , growth parameter b — d + c and activity c are related by a 
'subduality', i.e., a duality formula with a nonnegative error term. The next proposition generalizes 
this to branco-processes with positive annihilation rate. 

Proposition 17 (Subduality with a branching process) Let X be the {q,a,b,c,d)-branco- 
process and let y be the super random walk with underlying motion q^ , growth parameter 2a+b — d+c 
and activity 2a + c. Then 

E^[g-(0,Xt)j >iE0^g-(3^t,.T)j {xe£^{A), ^e[0,oof, \(t>\<oo). (2.73) 

Proof We first prove the statement if Ixl < 00 and \(f)\ < cxo. This goes exactly in the same way as 
in the proof of [AS05| Prop. 23]. Let H denote the generator oiy, defined in jAS05| formula (4.14)], 
let G be the generator in ()1.2|) . and let ^ be the duality function *(x, (f>) := e~^'^'^^ Then one has 

G^i; cp){x) - H^>{x, ■){^) = { ^q(z, j>(*)(e*«-*(^) - 1 - (0(*) - ^(j'))) 

+aY,x(i){x{i) - l)(e2*W - 1 - 2^{i)) + 6^x(z)(e-*W - 1 -t- </>(*)) (2.74) 

i i 

+c^a;(^)(2;(^) - l)(e"^W _ 1 _ 0(i)) + d^ x(i)(e"*« - 1 - 0(^))}e-<'^'^) > 0. 

i i 

This is just |AS051 formula (4.19)], where the extra terms with the prefactor a obtain their e^"^*-*-* — 1 
part from the generator G and the remaining —2(f){i) from H. Using Proposition [5] to guarantee 
integrability we may apply |AS09bl Corollary 2] to deduce (|2.73p . 

To generalize (j2.73p to x G '£'7 (A) and (p G [0, 00) supported on a finite set, we choose A/'(A) 3 
Xn^ X and let X" denote the (g, a, 6, c, (i)-branco-process started in X^ — x„. Then Corollary [TT] 
implies that the X" can be coupled such that X^{i) -^ Xt{i) a.s. for each i e A. It follows that 
g-{'l>,x^) _^ g-(0,Xt) g^ g^ ^^^ g-(3^t,a:„> ^ g-(i't,x„) g^ g^^ gQ taking the limit in ((^751) we obtain the 

statement for x e ^7 (A) and (p finitely supported. To generalize this to \4>\ < cx) we choose (pn^ (p 
and let 3^" denote the super random walk started in y^ — (pn- Then it is well-known (and can 
be proved in the same way as [AS051 Lemma 22]) that the 3^" can be coupled in such a way that 
y" t yt for each t > 0. Therefore, taking the monotone limit in (|2.73p our claim follows. I 
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2.6 The process started at infinity 

In view of what follows, we recall the following projective limit theorem. Let E and {Ei)i^fi be 
Polish spaces. Assume that ni : E ^- Ei are continuous surjective maps that separate points, 
i.e., for all x,y d E with x ^ y, there exists an i G N with ni^x) ^ T^iiy)- For each i < j, let 
TTij : Ej — > Ei be continuous maps satisfying tt^j o ttj = ttj . Assume moreover that for each sequence 
{xi)i,zfi with Xi € Ei [i €i N) that is consistent in the sense that ■nij{xj) — Xi {i < j), there exists an 
y € E such that Ttiiy) — Xi (« S N). Let {pi)i^fi be probability measures on the Ei's, respectively 
(equipped with their Borel-cr- fields) , that are consistent in the sense that /i^ — jij ott~j for all i < j. 
Then there exists a unique probability measure ii on E such that fi o tt^ — /ij for all i G N. 

This may be proved by invoking Kolmogorov's extension theorem to construct a probability 
measure ^' on the product space Yli Ei whose marginals are the Hi and that is moreover concentrated 
on the set E' C Yii ^i consisting of all {xi)i^fq satisfying Trij{xj) = Xi for all i < j. Now 7r(j/) := 
(7ri(?/))igN defines a bijection tt : i? — >■ £", so there exists a unique measure fi on the cr-algebra 
generated by the {TTi{x))i^m whose image under tt equals fi' . By |Sch731 Lemma IL18], this a- 
algebra coincides with the Borel-cr-algebra on E. 

Proof of Theorem[4]In the case without annihilation, parts (a)-(e) were proved in |AS051 Thm 2]. 
The proof there made essential use of monotonicity, which is not available in case a > 0. Instead of 
trying to adapt these arguments, replacing monotone convergence by some other form of convergence 
wherever necessary, we will make use of Corollary |3l which will simplify our life considerably. 

In view of this, set a := a/{a + c) and let X be the {q, 0, {l + a)b,a + c, a5 + d)-branco- process 
started at infinity, as defined in |AS051 Thm 2]. Fix e > and let {Xf)t>e be a {q, a, b, c, ci)-branco- 

process started at time e in A"| — Thin^_(Ar^ ). It has been proved in |AS05[ Thm 2] that 

X|°°^ G f^(A) for all t > a.s., hence Thin_L(A^°°') G f^(A) and {X[)t>e is well-defined by 
Proposition [9] By Corollary [3l 

¥[XI G •] =P[Thin_^(A|°°^) G •] (i > e). (2.75) 

In particular, this implies that if we construct two processes X^, X^ for two values < e < e', then 
these are consistent in the sense that {X^)t>£' is equally distributed with (A"f )t>e'- By applying 
the projective limit theorem sketched above, using the spaces of componentwise cadlag functions 
from (e, oo) to N^, we may construct a process {Xf')t>a such that X^ is equally distributed 
with Thin^^(A"^ ) for all e > and (A^ )i>o evolves as a (g, a, 6, c, (i)-branco-process. Let 

X'' denote the (g, r, s, ?Ti)-resem process with r,s,m as in (|1.8p . Then X'^ is dual to both the 
(q, a, &, c, (i)-branco-process (with parameter a — a/ [a + c) in the duality function) and to the 
(g, 0, (1 + a)&, a + c,ah-\- d)-branco-process (with duality function *(x, </>) = (1 — (t))''^)- We have 

\-(°o) Thin 1 ('T*-°°^ -E^(o°) 

E[(l-(l + a)0)^* ]=E[(l-(l + a)</>) tt^^ * ^] =^[{l - <t>)^t ] ^^.TG) 

= f^[x}=Q] (i>0, 0G[O,1]^, 101 <oo), 

where the last equality follows from |AS05| formula (5.5)] and we assume |0| < oo to make sure 
the infinite products are well-defined. It has been shown in |AS05| Thm 2 (d)] that the law of 

A( converges as i ^- oo to an invariant law of the (g, 0, (1 -I- a)h, a + c, ab + (i)-branco-process. 
By thinning, it follows that the law of X\°° converges as t — >■ oo to an invariant law V of the 
(g, a, 6, c, (i)-branco-process. Taking the limit f — ^ oo in (|2.76p we arrive at (|1.22p . Setting 

r ■.= {l + a)b+[a + c)-[ab + d) = a + b + c-d, (2.77) 

we obtain from |AS05[ Thm 2 (b)] and the fact that A(°°) is a 1/(1 + a)-thinning of A^""-*, that 

E[A^)(^)] < I ^(^^*^ ;['^ ^ JJ' (, e A), (2.78) 

L T+^ I^I+Sjt It r - U 
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which by the fact that 1/(1 + a) = {a + c)/(2a + c) yields ([TTH]). Formula pT^ is a simple 

consequence of the way we have defined X^°°'> as a thinning of X . This completes the proof of 
parts (a), (b), and (d)-(f) of the theorem. 

To prove also part (c), by formula (|2.76p and duality, it suffices to show that for each t > 

E[(l-(l+«)0)^* ]=E^[{l~{l+a)xlf ]^¥[X^^0] (0 e [0, 1]^, |0| < ^). (2.79) 

By Lemma [21] (i) below, X^ (i) < 1 a.s. for all i G A, hence a.s. on the event A"/ 7^ there exists 
some i e A such that < Xi (i) < 1. It follows that |1 — (1 + a)<Y/|^ ^> as n — >■ 00 a.s. on the 
event that A^ 7^ 0, hence (|2.79p follows from bounded pointwise convergence. I 

Remark Let X^") be as in Theorem|3](c). Then, using Proposition [T71 copying the proof of [AS05[ 
Thm 2 (b)], we obtain the uniform estimate 

if r' ^ 0, 



E[X^t'^\i)] < \ l^« + cjii-e ) ^.^^^^ ^2.80) 

77^ r if r' = 

{2a + c)t 

where r' := 2a + b + c — d. It is easy to see that this estimate is always worse than the estimate 
(jl.lSp that we obtained with the help of thinning (Corollary ^ . 

2.7 Particles everywhere 

The aim of this section is to prove Lemma ^U\ below, which, roughly speaking, says that if we 
start a (g, a, b, c, (i)-branco-process in a nontrivial spatially homogeneous initial law, then for each 
i > 0, if we look at sufficiently many sites, then we are sure to find a particle somewhere. For zero 
annihilation rate, this has been proved in |AS05[ Lemma 6] . Results of this type are well-known, see 
e.g. the proof of |Lig85[ Thm III.5.18]. It seems the main idea of the proof, and in particular the use 
of Holder's inequality in (|2.88|) below or in |Lig85[ (III. 5. 30)] goes back to Harris [Har76| . Another 
essential ingredient of the proof is some form of almost independence for events that are sufficiently 
far apart. For systems where the number of particle per site is bounded from above, such asymptotic 
independence follows from Lig85[ Thm 1.4.6], but for branco-processes, the uniform estimate given 



there is not available. In |AS05j . we solved this problem by using monotonicity, which is also not 
available in the presence of annihilation. Instead, we will base our proof on the covariance estimate 
from Proposition [13] above. 

Lemma 18 (Particles at the origin) Let G be a transitive subgroup 0/ Aut(A) and let fi be a 

G-homogeneous probability measure on £^(A). Assume that b > 0. Then, for a.e. x w.r.t. ^, the 
(g, a, 6, c, d)-branco-process started in Xq = x satisfies 

P^[Xt(0) >0] >0 (t>0). (2.81) 

Proof Although the statement is intuitively obvious, some work is needed to make this rigorous. If 
a ~ 0, then by monotonicity (see LemmalU which extends to infinite initial states by Corollarv fTT]) . 
it suffices to prove that for a.e. x w.r.t. /j,, there exists some i S A with x{i) > such that there is 
a positive probability that a random walk with jump rates g, started in z, is at time t in the origin. 
Since we are only assuming a weak form of irreducibility (see (II. ip (ii)), this is not entirely obvious, 
but it is nevertheless true as has been proved in [AS05[ Lemma 31]. 

If a > 0, then, to avoid problems stemming from the non-monotonicity of X, we use duality. Let 
a, r, s, m be as in (|1.8p and observe that ttt, > by our assumptions that a, 5 > 0. Define 5q S [0, 1]"^ 
by 5[){i) :— 1/^=0} • Then, by duality (Proposition[T|), letting X denote the {q\ r, s, m)-resem-process 
started in Xq = 6q, we have 

E"[(l-(l + a))^*(0)] =E*'«[(l-(l + a)A'/)"], (2.82) 
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and our claim will follow once we show that for all t > 0, this quantity is strictly less than one 
for a.e. x w.r.t. /i. Thus, it suffices to show that P''°[0 < Xt{i) < 1] > for some i 6 A such that 
x{i) > 0. By the fact that m > and Lemma \T\\ (i) below, this can be relaxed to showing that 
V^°[Xt{i) > 0] > for some i E A such that x{i) > 0. Letting X denote the (g, 0, s, r, m)-branco- 
process, using duality again (this time with a = 0), it suffices to show that 

1>E-^«[{1-Xlf] =E^[0^*('^)] =P^[lt(0)=0]. (2.83) 

Thus, the statement for systems with annihilation rate a > follows from the statement for systems 
with a = 0. I 

Lemma 19 (Finiteness of moments) Let X be a {q, a, b, c, d)-branco-process started in an arbi- 
trary initial law on £j{A). Assume that (A,g) is homogeneous and that a + c > 0. Then 

E[Xi(i)"] <cx) (to> 1, ie A, i>0). (2.84) 

Proof By Lemma [7] and Corollary [TTl for each i > we can couple a (g, a, 6, c, (i)-branco-process 
X started in an arbitrary initial law on £-y{A) to the {q,0,b,a + c, d)-branco- process X' started 
in the same initial law, in such a way that Xt < X^ a.s. In view of this, it suffices to prove the 
statement for the system X' with zero annihilation rate and annihilation rate c' := a + c. Let 
X'(") be the (g, 0, fe, c',d)-branco-process started in X'^"'>{i) = X^(i) V n (i G A). Then, by |AS05[ 

Theorem 2 (c)], for each i > the process X^' can be coupled to the process started at infinity, 
denoted by X^°°\ in such a way that ATj -^ t Xl°°' a.s. In view of this, it suffices to prove that for 
the process without annihiation started at infinity 

E[X^^°°\i)"'] < oo (to > 1, i G A, i > 0). (2.85) 

Let A"(") denote the (q, 0, &, c',d)-branco-process started in the constant initial state Xq (i) = n 
{i e A). Again by |AS05[ Theorem 2 (c)], it suffices to find upper bounds on E[X^/'\i)"'] that 
are uniform in n. Such upper bounds have been derived in [DDLQOl Lemma (2.13)] for branching- 
coalescing particle systems on Z'^ with more general branching mechanisms than considered in the 
present paper. In particular, their result includes (g, 0, &, c',(i)-branco-processes on Z*^ with c' > 0. 
Their arguments are not restricted to 'E'^ and apply more generally to underlying lattices A and 
jump kernels q as considered in the present paper, as long as (A, q) is homogeneous. I 

Remark It seems likely that the assumption in Lemma[T9]that (A, q) is homogeneous is not needed. 
The proof of |DDL901 Lemma (2.13)], which we apply here, uses translation invariance in an essential 
way, however. Since we do not need Lemma[T2]in the inhomogeneous case, we will be satisfied with 
the present statement. It does not seem easy to adapt the proof of formula (|1.18p (which holds 
without a homogeneity assumption) to obtain estimates for higher moments. 

Lemma 20 (Systems with particles everywhere) Assume that (A,q) is infinite and homoge- 
neous, G is a transitive subgroup o/ Aut(A, q), and a-\-c > 0, b > 0. Let X be a (g, a, 6, c, d)-branco- 
process started in a G-homogeneous nontrivial initial law on £y{A). Then, for any t > 0, 

lim P[ThinA (Xf) = 0] = (2.86) 

n— >oo 

for all (f)n G [0, 1] satisfying \(t>n\ ^ oo. 

Proof By Lemma [T9l restarting the process at some small positive time if necessary, we can without 
loss of generality assume that E[Xo(0)^] < cxd. Set 7r„ := 4'n/\4>n\ and let P^ denote the law of the 
process started in a deterministic initial state x. Then, for each r < oo and i > 0, we can choose n 
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sufficiently large such that r < \4>n\- Then a r7r„-thinning is stochastically less than a (/)„-thinning 
and therefore 

P"[Thin^„(XO - 0] < P"[Thin,,„(XO = 0] 

where in the last step we have applied Jensen's inequality to the concave function z i— )■ z^^"^"^' . For 
the process started in a nontrivial homogeneous initial law, we obtain, using Holder's inequality, 
for all n sufficiently large such that r < |0„|, 

P[Thin^„(Xt) - 0] = fp[Xo e da;]P"[Thin0„(Xt) = 0] 

< n ( fnXo e dx]E-[e-^*(*)]y ^%E[i?„(Xo)] (2.88) 

= n ( /p[^oedx]E-[e-^*(0)]y"^'VEK(Xo)] 

= /"p[Xo e dx]E" [e-^tW]" + E[i?„(Xo)], 

where we have used spatial homogeneity in the last step but one. 

By CoroUarvfTSl the quantity Rn{x) defined in (|2.87p can be estimated as 

\Rn{x)\<\r'' J2 (5I-'^W-^*(*;^'0 + I]a:«xO-)it(z,j;fc,0)^„(fc)7r„(0. (2.89) 

k,i * ''^ 

kjtl 

It follows that 

E[|i?„(Xo(0))|] <y Y. {Y.^[M^)]Kt{^;kJ) + YE[Xo{i)Xo{J)]Lt{^,J■,k,l)y,,{k)7^r,{l) 
k,i * '^3 

<y Yl (E[^o(0)]^7^t(*;fc,0 + E[^o(0)']^£t(*,i;fc,0)^n(fc)7r„(0 

k,l ■■' 

k=^l 

=:r2^C(fc,0^„(fc)7r„(/), 

k,l 

(2.90) 
where by definition C(fc, k) := and we have used that by Cauchy- Schwartz and translation invari- 
ance: 

|E[Xo(*)Xo(j)]| < E[Xo(*)2]i/2E[Xo(jf ]i/2 = E[Xo{On (2.91) 

We claim that 

y Cik,l)7Tn{k)TT„il) -^ 0. (2.92) 

k,l 

To see this, we observe that by (|2.34p . (|2.35p and our assumption that E[Xo(0)^] < oo, 

C{gk, gl) = C(fc, /) {g £ G) and ^ C{k, 0) < cx). (2.93) 
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Since G is transitive, for each Z S A we can choose some gi (z G such that gil = 0. In view of this, 
(|2.93p shows in particular that for each e > 0, the quantity 

|{fc e A : G{k,l) > e}\ = \{gik e A : C{gik,0) > e}\ = \{j G A : C(j, 0) > £}| =: K, (2.94) 

does not depend on / e A and is finite. It fohows that 

^7r„(0^C(fc,ZK„(fc)<^7r„(0( J2 ''n{k)+ ^ 7r„(fc)) < if,/|(/)„| + e. (2.95) 

I k I k:C(k,l)>e k:C{k,l)<e 

Since |(/)„| -^ oo and e > is arbitrary, this proves (|2.92p . By (|2.88p and (|2.90p . we conclude that 
for each r < oo, 

linisupP[Thin0„(Xt) = 0] < f P[Xo G dx]E"= [e^^* Wj''^ (2.96) 

Letting r — > cx), using 6 > and Lemma [T51 we arrive at (|2.86p . I 

2.8 Long-time limit law 

In this section, we prove Theorem [5] We first need some preparatory results. 

Lemma 21 (Not exactly one) Let X be a {q, r, s, m)-resem process started in a finite initial state 
(j> G [0, 1] , |0| < oo. Assume that (A, g) is infinite and homogeneous and that m > 0. Then 

(i) P'''[Xt{i) = 1] = for each t > 0, i G A. 

(ii) P'*[0 < \Xt A (I - Xt)\ < K] ^ as t^ oo for all K <oo. 

Proof Let X~^ and X~ satisfy ^Yq^ = Xq = Xq — (j) and be given, for times i > 0, by the solutions 
to the stochastic differential equations 

dX+it) = {\q\ +s)(l- X+{t))dt - mX+{t)dt + J2rX+{i){l ~ X+ {i))dBt{i) 

I (2.97) 

dX^{i)^^{\q\^m)X^(i)dt+^2rX^{i){l-X^-(i))dBt{i) (t > 0, i G A), 

where \q\ := '}2ijq{i,j), which does not depend on j G A by the transitivity of Aut(A,q), and 
(i3(i))igA is the same collection of independent Brownian motions as those driving X. By the 
arguments used in the proof of [AS05[ Lemma 18], solutions of (|2.97p are pathwise unique and 
satisfy 

Xf- <Xt< X+ {t > 0) a.s. (2.98) 

Moreover, since (|2.97p contains no interaction terms, the [0, l]^-valued processes {Xf {i),Xf^{i))t>o 
are independent for different values of i G A. Since X'^{i) is a one-dimensional diffusion with (by 
grace of the fact that tti > 0) the drift on the boundary point 1 pointing inwards, it can be proved 
by standard methods that 

f[X+{i) = 1] = (i > 0, i G A). (2.99) 

We defer a precise proof of this fact to Lemma [23] in the appendix. Together with (|2.98p , formula 
(|2.99|) proves part (i) of the lemma. 

To prove also part (ii), we observe that 

E[A'r(i)]-e-(l«l+™)V(i) (t>0, iGA). (2.100) 

With a bit of work, it is possible to show that there exists a io > and function (0, tg] 3 1 1-^ ct > 
such that 

E[Xf{i) A (1 - X+{i))] > ct(t>{i) {0<t< to). (2.101) 
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A precise proof of this fact can be found in Lemma [25] of the appendix. We note that for any [0, 1]- 
valued random variable Z, one has Var(Z) = E[{Z -E[Z])'^] < E[\Z -E[Z]\] < E[Z + E[Z]] = 2E[Z]. 
Applying this to Z = -%"("(«) A (1 - A;+(z)), using (|2.100p . we see that 

Yai{X-{i) A (1 - X+{i))) < 2e-(l«l +"^)V(i) (i > 0, i€ A). (2.102) 

Now (|2.10ip implies E[|A'(^ A(l-A:'f+)|] > ct\4>\, while by (|2.102l) and the independence of coordinates 

i e A, 

Var(|A'j-A(l-A'+)|) <2e"(l'?l+™)Vl (0 < t < io)- (2.103) 

Since Xt{i) A (1 ~ Xt{i)) > Xf{i) A (1 - A;+(i)), by Chebyshev, it follows that 

P^[\XtA{l-Xt)\<h,Ct\<P\] < ^ l*^' (0<t<to), (2.104) 

which tends to zero for \4>\ — > oo. By |AS05[ Lemma 5], 

P'l'\0<\Xt\<K] — >0 {K<oo). (2.105) 

It follows that we can choose Lj — > oo slow enough such that 

P'l'\0< \Xt\ < Lt] — > 0. (2.106) 

By (|2.104p . we conclude that 

limsupP"*[0 < \Xt A{l-Xt)\< K] 



(2.107) 



< limsupP*[0 < \Xt A{l-Xt)\<K\0< \Xt-to\ < it-to]P^[0 < \Xt^t„\ < Lt-t, 
+ limsupP'^[0 < \Xt A (1 - Xt)\ <K\\Xt-t„\ > Lt-to]P^[\Xt-to\ > Lt-t,] 

< limsupP'^[|A't A{l-Xt)\< \ct,Lt_t, \ \Xt-to\ > Lt-t,] 

< limsup ^ ^ „ = 0. 

■ 

Remark It seems likely that the condition m > in Lemma [^ is not necessary, at least for part (i). 
Indeed, it seems likely that ((7,r, ,s,m)-resem-processes have the 'noncompact support property' 

F[Xt{i) > 0, Xtij) - 0] = (i > 0, i,j G A, q{i,j) > 0), (2.108) 

similar to what is known for super random walks |EP91I . Since proving (|2.108p is quite involved 
and we don't know a reference, we will be satisfied with proving Lemma [21] only for m > 0, which 
is sufficient for our purposes. 

Lemma 22 (Systems with particles everywhere) Assume that (A, q) is infinite and homo- 
geneous and that G is a transitive subgroup of Aut(A, g) and a + c > 0, b > 0. Let X be the 
(g, a, 6, c, d)-branco process started in a G-homogeneous nontrivial initial law £(^0). Then, for any 
t > and < a < 1 and for any e > 0, there exists a K < 00 such that 

|(/)| < 00 and I0A (1-0)1 > K implies E[\1 - (l + a)0|^*] < e. (2.109) 

Proof We start by proving that if 4>n S [0, 1]^ satisfy |(/)„| < 00 and |0„ A (1 — 0„)| — > 00, then 

lim E[|l-(l + a)0„|-^*] =0. (2.110) 
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Set i/jn '■— 4>n A (1 ~ (f)n)- Then, for each i e A, we have and ipn{i) < ^ — 4>n{i) < 2 — (1 + «)(/)„(«) 
and ipnii) < '/'«(«) < (1 + 0!.)(j)n{i), from which we see that 

Mi) - 1 < 1 - (1 + a)c^n{i) < 1 - Mi): (2.111) 

or, in other words, |1 — (1 + a)(/)„(i)| < 1 — ipn- It follows that 

|E[l-(l + a)0^*]| < E[|l-(1 + «)(/.„ 1^'] <E[(l-7/^„)^'] =P[Thinv,„(Xt)-0], (2.112) 

which tends to zero by Lemma BUI and our assumption that {ipnl -^ c»- 

Now imagine that the lemma does not hold. Then there exists some e > such that for all 
n > 1 we can choose (/)„ with |0„| < cxd and |0„ A (1 — 0„)| > n such that E[|l — (1 + a)0„|'^'] > e. 
Since this contradicts (|2.110p . we conclude that the lemma must hold. I 

Proof of Theorem [5] For a = the statement has been proved in |AS05| Thm 4 (a)], so without 
loss of generality we may assume that a > 0. By Theorem |4] (e) , it suffices to show that 

E[(l-(l + a)0)^*l -^ P"^ [3i > such that A-/ = 0] {(pe[0,lf, |(/)| < oo), (2.113) 

where a :— a/ {a + c) and X'^ denotes the (g^, a + c, (1 + a)5, ab + d))-resem-process started in (f). 
By duality (Proposition [1]) , for each t > 1, 

E[(l - (1 + a)<i>f'] = E[(l - (1 + a)xl,f'] , (2.114) 

where X'^ is independent of X and started in Xq = (j). For each K < oo, we may write 

E[{l^{l + a)xl,f']^F[xl,^0] 

+E[(1 - (1 + a)Xl,)^' I < \Xl, A (1 - Xl,)\ < K]¥[0 < \xl, A (1 - xl,)\ < K] 

+E[(1 - {l + a)xi,)^' I K < \xi, A (1 - XI,)\]F[K < \xl, A (1 - A-.l^)!]. 

(2.115) 
Here the first term converges, as i — ?> oo, to P"^ [3t > such that X^ = O] . Note that a > by our 
assumption that a > 0. Assume for the moment that also 6 > 0. Then Lemma [H] (ii) tells us that 
the second term on the right-hand side of (I2.115P tends to zero. By Lemma [22l for each e > we 
can choose K large enough such that the third term is bounded in absolute value by s. Putting 
these things together, we arrive at (|2.113p . 

If fe = 0, then Lemma [?11 (ii) is not available, but in this case \Xf\ is a supermartingale, hence 
|AS051 Lemma 5] tells us that P*^ [Eli > such that A"/ = O] =1, and the proof proceeds as above.! 

A Some facts about coupled Wright-Fisher diffusions 

The aim of this appendix is to prove two simple facts about (coupled) Wright-Fisher diffusions. 
In particular, applying Lemmas [23l and [26l to X — X^{i), Y = X^{i), a = \q\ + s, b = rn and 
c = I?! + m yields formulas (|2.99p and (|2.10ip . respectively. 

For a, 6 > and r > 0, let X denote the pathwise unique (by |YW71) ) [0, l]-valued solution to 
the stochastic differential equation 



dXt = a(l - Xt)dt - bXtdt + y/2rX{l - X)dBt, (A.l) 

where B is standard Brownian motion. 
Lemma 23 (No mass on boundary) If b > 0, then 

P[Xt = 1] = it> 0), (A.2) 

regardless of the initial law. 
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Proof If a, 6 > 0, then it is well known that X has a transition density (see Propositions 3 and 4 in 
[Palll] along with the discussion on page 1183 or [ Gri79b[ IGri79a| V Consequently ¥[Xt = 1] = 
and hence the result follows. If a = but b > 0, then by standard comparision results (see |Bas98[ 
Thm. 6.2] or IAS05[ Lemma 18]), if Zq = Xa and Z solves the SDE (fOt with a = 6/2 and b replaced 
by b/2, relative to the same Brownian motion, then Xt < Zt and hence F[Xt = 1] < P[^t = 1] = 
for all i > 0. I 

Lemma 24 (Moment dual) Let K = {Kt)t>o be a Markov process with state space N U {oo}, 
where oo is a trap, and K jumps from states k eN as 



fc I— > A; — 1 with rate ak + rk{k — 1), 

fc I— ). oo with rate bk. 



(A.3) 



Ther 



E'^iX^] = E''[x'^'] {t > 0, xe [0, 1], k e N), (A.4) 

where x'^ := 1 and x°° := for all x G [0, 1]. 

Proof Let 

Gfix) := [a{l x)l - bxl + rx{l - x)^] fix), 

Gf{k) := [ak + rk{k - 1)] {/(fc - 1) - /(fc)} + fefc{/(oo) - f{k)] 

be the generators of the processes X and K, respectively, and let il^{x, k) := x'' be the duality 
function. Then 

gip{ ■ , k){x) = ak{x''-^ ~ x'') - bkx'' + rk{k - l){x''-^ - x'') ^ GV'(x, , • ){k), (A.6) 

where the term with k{k — 1) is zero for k — 1 and both sides of the equation are zero for fc = 0. 
The claim now follows from [AS05[ Thm 7] and [ASOQb] and the fact that the expression in (|A.6p 
is bounded uniformly in x and k, which guarantees the required integrability. 

Although this is not needed for the proof, this duality may be understood as follows. We can 
view Xt as the frequency of type-one organisms in a large population where pairs of organisms are 
resampled with rate 2r and organisms mutate to type 1 and 0, respectively, with rates a and b. Then 
E[X('^] is the probability that k organisms, sampled from the population at time t, are all of type 
one. We can view Kt as the ancestors of these organism at time zero, where we neglect organisms 
that due to mutation are sure to be of type one while on the other hand the state Kt = oo signifies 
that due to a mutatation event, at least one of these ancestors is of type zero. I 

Now let X be as in (|A.1|) . let c > 0, and let Y be given by the pathwise unique solution to the 
stochastic differential equation 



dYt = -cYtdt + y/2rY{l ~ Y)dBt, (A.7) 

driven by the same Brownian motion as X. 

Lemma 25 (Feller property) Let {X, Y) be given by the pathwise unique solutions of HA . 1]} and 
IjA.T^ I, and let Kt{{x,y), ■) :— F^^'y^[{Xt,Yt) G •] denote the transition probabilities of {X,Y). 
Then the map {t, x, y) M- Kt{{x, y), ■ ) from [0, oo) x [0, 1] into the probability measures on [0, 1]^ is 
continuous w.r.t. weak convergence of probability measures. 

Proof It follows from well-known results |EK861 Corollary 5.3.4 and Theorem 5.3.6] that pathwise 
uniqueness for a stochastic differential equation implies uniqueness of solutions to the martingale 
problem for the associated differential operator, which is in our case given by 



A := a(l - x)^ - bx^ + rx{l - x)^ - cy^ + ry{l - y)|^ + 2r^ xil ^ x)y{l ~ y)^, (A.8) 
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(A.11) 



with domain C^[0, 1]^. Now if {X", Y") are solutions to this martingale problem with deterministic 
initial states {Xq,Yq) = {xn,yn) converging to some limit {x,y) S [0, 1]^, and {X,Y) denotes the 
process started in {x,y), then |EK86[ Lemma 4.5.1 and Remark 4.5.2] imply that 

P[(xr,yr)*>o e •] ^ ¥[{Xt,Yt)t>o e •], (A.9) 

where => denotes weak convergence of probability laws on the space C[o_i]2[0, oo) of continuous 
functions from [0,oo) into [0,1]^, equipped with the topology of locally uniform convergence. In 
particular, this implies the stated continuity of the transition probabilities. I 

Lemma 26 (Linear estimate) Assume that b > 0. Then there exists a to > and function 
(0, io] 3 1 1-^ Xt > such that the process started in {Xq, Yq) — {z, z) satisfies 

E'^^'^^[YtA{l-Xt)]>XtZ (0<i<io, 0<z< 1). (A.IO) 

Proof We estimate 

E(^'^) [Ft A (1 - X,)] >E(^^^)[Ft(l-XO] >iE(^''^[5^tl{jf,<i}] 
- l{E^[Yt]-W[l^x,>i}]) > KM -2IE^[^']' 
where the last step we have used that lr3.>ij < 4a;^. By Lemma [24l 

E'[Yt] = e-^^z {t>0, z e [0, 1]), (A.12) 

while by the same lemma 

E^\X?] ^E^fz-^'l <F^\Kt < l]z + ¥^\Kt ^ 2\z^ 

L I J (A 13) 

^l^l_^-2ia+r)t^^^^^2ia+b+r)t^2^2{a + r)tz + z\ 

Combining this with (jA.ll[) yields 

E(^.^) [Yt A (1 - Xt)] > {^e-"' ~ 4(a + r)t - 2z)z. (A.14) 

Choosing ig > and zq > small enough, we find that 

E(^'^) [Yf A (1 - Xf)] > iz (0<t <io, 0<z< zo)- (A.15) 

To extend this to all z £ [0,1], at the cost of assuming that t > and replacing the constant 
1/4 by a possibly worse, time-dependent constant At, we observe that by Lemma [^ the function 
[0, l]3 z^ E(^'^) [Yt A (1 - Xt)] is continuous. Since by Lemma [Hand (|XT2l) . 

E(^'^)[rt A(l-Xt)] >0 (t>0, ze (0,1]), (A.16) 

using continuity, we may estimate E*^^'^^ [Yj A (1 — Xt)\ uniformly from below on [zo,l], which 
together with (|A.15|) yields (|A.10|) . I 
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